Volume behaviour of quark condensate, pion mass and decay constant from 

Dyson- Schwinger equations 
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We solve the coupled system of Dyson- Schwinger and Bethe-Salpeter equations for the quark 
propagator and the pion Bethe-Salpeter amplitude on a finite volume. To this end we use a trun- 
cation scheme that includes pion cloud effects in the quark propagator and light mesons. We study 
volume effects in the quark condensate, the pion mass and the pion decay constant and compare 
to corresponding results in other approaches. In general we find large effects for volumes below 
V = (1.8 fin) 4 . 
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I. INTRODUCTION 

Dynamical chiral symmetry breaking is one of 
the fundamental properties of QCD. Its pattern de- 
termines the experimentally observable spectrum 
of light hadrons as well as the details of the under- 
lying interaction between quarks and gluons. It 
therefore plays a central role in our understanding 
of QCD. Despite great efforts in the last decades 
not all details of dynamical chiral symmetry break- 
ing are as well understood as one could wish. 

In this respect, the chiral properties of QCD at 
finite volume have attracted considerable attention 
over the years. Strictly speaking, dynamical chi- 
ral symmetry breaking and the associated forma- 
tion of Goldstone bosons is restricted to the infinite 
volume limit. On any given finite volume the cor- 
responding order parameter, the quark condensate 
(qq), goes to zero in the limit of vanishing bare 
quark mass m. Nevertheless, one can extract the 
properties of the infinite volume theory from the 
formulation in a box. For large enough volumes 
and/or quark masses the box effects are small and 
the finite volume theory closely resembles its infi- 
nite volume limit. This happens if the condition 



to V (qq) > f 



(1) 



is satisfied and the eigenvalues of the Dirac op- 
erator are almost dense [l|. On the other hand, if 
to V (qq) <C 1 dynamical chiral symmetry breaking 
is lost. 

Another important condition is the relation of 
the four-volume V with a typical hadronic scale 
Ah and the pseudo-Goldstone mass M„. If the 
inequality 

k <<¥<< ^ (2) 

is satisfied, the theory is in a critical region where 
long range correlations are at work. Here the 



QCD partition function depends on the quark mass 
to and the volume only through the combination 
fi := mVY,, where S is the chiral condensate in 
the infinite volume limit. In this critical region ex- 
act, analytic and universal scaling laws are known 
from chiral random matrix theory see e.g. 

[j| for a recent review. Certainly these are of great 
importance for approaches like lattice QCD which 
are inherently limited to finite volumes and lattice 
spacings. 

The properties of QCD at finite volumes have 
been studied by several approaches. One, of 
course, is lattice QCD, see @-[lll 

and Refs. 

therein, but important insights have also been ob- 
tained from effective theories such as chiral pertur- 
bation theory fl2l - [l4j ]. or the quark- meson model 
In this work we follow a complemen- 
tary approach. We study volume effects in the 
quark condensate, the pion mass and the pion 
decay constant from the framework of Dyson- 
Schwinger equations (DSEs) of the quark propa- 
gator and the corresponding Bethe-Salpeter equa- 
tions (BSEs) for pseudoscalar bound state ampli- 
tudes. 

One of the advantages of this approach is that 
volume effects can be studied continuously from 
very small to very large volumes (corresponding 
studies for meson observables using chiral pertur- 
bation theory for example have to distinguish be- 
tween two different regions of chiral counting). The 
implementation of mixed boundary conditions in 
the spatial and time directions is possible without 
great effort. Furthermore one has direct access to 
the infinite volume and the continuum limit with- 
out the need to perform any extrapolations. We 
are thus in a position to study chiral symmetry 
restoration at small volumes together with effects 
at large and infinite volumes in the same frame- 
work. 



FIG. 1: The Dyson-Schwinger equation for the quark propagator. Filled circles denote dressed propagators 
whereas the empty circle stands for the dressed quark-gluon vertex. 



On the other hand we need to work with a trun- 
cation scheme in order to close our system of equa- 
tions. This scheme is systematic in the following 
sense: in a first step we neglect the volume depen- 
dence of the three-point functions. For the quark- 
gluon vertex we work with a volume independent 
ansatz that has proven useful in calculations of 
pion cloud effects in the masses of light mesons 
HH EH • The volume dependence of the gluon and 
quark propagators arc then determined from their 
Dyson-Schwinger equations and serve as input into 
the volume dependent Bethe-Salpeter equation of 
the pion. This truncation scheme can be system- 
atically extended to also include volume effects of 
higher Green's functions by explicitly solving their 
corresponding DSEs. 

This work is organised as follows: In the next 
two subsections we summarise the properties of the 
quark Dyson-Schwinger equation and the Bethe- 
Salpeter equations on a torus. We then discuss de- 
tails of the numerical treatment of these equations 
in subsection III CI whereas in subsection III Dl we 
specify our approximation scheme for the quark- 
gluon vertex. Our numerical results for the quark 
condensate are presented in section ITlI Al We then 
proceed to discuss the volume dependence of the 
pion mass and decay constant in subsection 1111 Bl 
Starting with subsection IIII Cl we also take into ac- 
count finite volume effects in the gluon propagator. 
We conclude and summarise in section ITVl 



II. QUARKS AND PIONS AT FINITE 
VOLUME 



A. The quark Dyson-Schwinger equation on 
a torus 

In Euclidean momentum space, the renormaliscd 
dressed gluon and quark propagators in the Lan- 
dau gauge are given by 



S(p) = 



p fl p„\ Z(p 2 ) 



1 

itfA{p 2 ) + B{p 2 ) 



(3) 
(4) 



Here the gluon dressing function Z(p 2 ) and the 
quark wave function Zf(p 2 ) = A^ 1 {p 2 ) also de- 



pend on the renormalisation point, whereas the 
quark mass function M(p 2 ) = B{p 2 )/A(p 2 ) is a 
renormalisation group invariant. These propaga- 
tors can be calculated from their Dyson-Schwinger 
equations (DSEs). Most important for this work is 
the DSE for the quark propagator which is shown 
diagrammatically in Fig. [TJ 

On a compact manifold, the gluon and quark 
fields have to obey appropriate boundary condi- 
tions in the time direction. These have to be pe- 
riodic for the gluon fields and antiperiodic for the 
quarks. It is convenient, though not necessary, to 
choose the same conditions in the spatial direc- 
tions. For the gluon propagator we will use pe- 
riodic boundary conditions in all four space-time 
directions. For the quark field we will employ an- 
tiperiodic boundary conditions in all four direc- 
tions. In principle it is also possible to implement 
periodic boundary conditions for the spatial direc- 
tions of the quarks; corresponding results will be 
presented elsewhere. 

We choose the box to be of equal length L 
in all spatial directions and a potentially longer 
value T in the time direction. The correspond- 
ing volume is denoted V = T L 3 . Together with 
the boundary conditions the finite volume leads 
to discretised values in momentum space. Thus 
all momentum integrals appearing in the Dyson- 
Schwinger equations are replaced by sums over 
Matsubara modes. For the ghost and gluon DSE 
the corresponding equations as well as their solu- 
tions have been discussed in detail in Rcfs. 
l22j . The corresponding solutions in the infinite 
volume/continuum limit are given and discussed 
in Ref. [23| . To make this paper self-contained we 
shortly repeat the general features of these solu- 
tions below, where we discuss the quark-gluon in- 
teraction. 

The quark-DSE on a torus has been discussed in 
Ref. [U, [25[ at zero and in Ref. [26| at finite tem- 
perature, where results for the quark propagator 
and condensate at real quark momenta have been 
given. In the present work the quark propagator 
also serves as input into the Bethe-Salpeter equa- 
tion for pseudoscalar mesons. As we will see this 
also necessitates the determination of the quark 
propagator for complex momenta. 

Restricting our system to a finite volume results 




FIG. 2: The discretised momentum space that will be 
used: the momenta for antiperiodic boundary condi- 
tions in all directions. The circles indicate that only 
'full' hyperspheres are taken into account, see text for 
details. This procedure leads to reduced cutoff effects. 



in the Matsubara modes 



2tt 

Pi = -f-rii 



Pi 



2ir 



1/2) 



(5) 



for periodic and antiperiodic boundary conditions 
in direction i with length Li and rij G Z. With dis- 
cretised momenta an integral in momentum space 
becomes 



d 4 q 



1 



N max 
E 

m ,712 ,™3 ,fi4 — — N 



(6) 



where N defines the maximal Matsubara mode, 
and therefore acts as an ultraviolet cutoff. Further- 
more N max = N for periodic and N max = N— 1 for 



antiperiodic boundary conditions. Similar to the 
infinite volume case where one uses hypersphcri- 
cal coordinates we rearrange the summation over 
Matsubara frequencies in the following convenient 
way 



E = EE 



(7) 



ni,n2,n^,n4 — — N 



where i counts hyperspheres containing all mo- 
menta with the same absolute value and m counts 
individual momenta on each sphere. In the follow- 
ing we frequently use the notation 



E : =EE; 



(8) 



to account for the sum over all momenta qj. 

Fig. [5] shows how the discretised momentum 
space looks in two dimensions, with the circles in- 
dicating the hyperspheres that are counted by i. 
The hyperspheres that are not fully occupied for a 
given value of N are not used, i.e. hyperspheres 
that receive additional vectors when going from ./V 
to N + 1 are discarded. This procedure leads to 
an approximate restoration of rotational symmetry 
for large momenta. Consequently cutoff artefacts 
in the ultraviolet momentum regions are drasti- 
cally reduced (HI 0, Hz| . 

The Dyson-Schwinger equation for the quark 
propagator on a torus is then given by 



S-\ Pi ) = Z 2 [S^i)}- 1 -C F ^j^Y.^ 5 feO r»(qj,Pi) - qj ) . (9) 

3 j 

I 



where the factor Cf =4/3 stems from the colour 
trace and we have introduced a reduced quark- 
gluon vertex T v (q,p), by defining TQ U (q,p) = 
ig^T„(q,p). Here p and q denote the momenta of 
the two quark legs, whereas k = p — q is the corre- 
sponding gluon momentum. The bare quark prop- 
agator is given by = ij-p+m(A' 2 ), where 
to(A 2 ) is the unrenormalised bare quark mass and 
A represents an ultraviolet cutoff. The wave func- 
tion renormalisation factor Z 2 is determined in the 
renormalisation process. The ghost renormalisa- 



tion factor Z% will be absorbed in our truncation 
of the quark-gluon vertex which we discuss below. 
The quark mass function M(p 2 ) and the wave func- 
tion Zf(p 2 ) can be extracted from Eq. © by suit- 
able projections in Dirac-space. 

B. The Bethe-Salpeter equation for the pion 
on a torus 

The homogeneous Bethe-Salpeter equation 
(BSE) for flavour non-singlet mesons on a torus 
can be written as 



Kp(^p) = m tY, K °**yfav> p ) W) rw fe; p )%7)] 7 a ( 10 ) 
3 



4 



where K is the Bcthc-Salpctcr kernel. The mo- 



menta g H 



£P and q~ = q - (1 - £)P of the 



quark constituents arc chosen such that the total 
meson momentum is given by P = q + — q~ . Here 
the momentum partitioning parameter £ = [0,1] 
reflects the arbitrariness in the relative momenta of 
the quark-antiquark pair. Since all obscrvables are 
independent of £ we can choose £ = \ without loss 
of generality. The flavour content of the meson is 
expressed through flavour matrices which are sup- 
pressed in Eq. (fTOj) . The Greek indices (a . . .) refer 
to colour and Dirac structure. The BSE is a para- 
metric eigenvalue equation with discrete solutions 
P 2 = —M 2 where M n is the mass of the bound- 
state. The lowest mass solution corresponds to the 
physical ground state. Since P 2 is negative, the 
momenta q ± are necessarily complex in Euclidean 
space and so the quark propagator functions must 
be evaluated with complex argument. This leads 
to technical issues that are dealt with in subsection 

ESI 

In general the pion Bethe-Salpeter amplitude 
can be decomposed into four different tensor struc- 
tures Pi.. 4 according to 

1"^; P) = 75 [Pi( Pl ; P) - iPF 2 ( Pl ; P) (11) 

-*/< {pi ■ P) F 3 ( Pl ; P) - P 4 (p 4 ; P) . 

The specific values for the momenta pi depend on 
our choice of boundary conditions for the quark 
fields. In the time direction the relative momenta 
Pi of the pion constituents have to be antfperiodic 
as can be seen from the exact expression |28| 



Fi( Pi ;P) = 



B( Pl 



(12) 



in the chiral limit, where B is the scalar quark 
dressing function and the pion decay constant. 
With our choice of antiperiodic boundary condi- 
tions for the quark fields in all four direction we 
obtain the same boundary conditions for the pi. 
The total momentum P of the pion is fixed by the 
condition — P 2 = AI 2 on the momentum shell of 
the pion. 

The (normalised) Bethe-Salpeter amplitudes of 
the pion can be used to determine the pion 
decay constant according to 



3Z 2 1 
M 2 L 3 T 



Y / tT[T^(q j ,-P)S(qf) 75 PS(qJ 



(13) 

on our torus. Here the trace is over Dirac matri- 
ces, and q ± = q ± P/2. As for the normalisation 
condition for the pion amplitude we use a prescrip- 
tion originally proposed by Nakanishi in Ref. [29( 
and put to use for momentum dependent kernels 
in 13011 . 



C. Solving the DSE for complex momenta 

As we have seen, in the pion BSE the quark 
propagators are required for complex momenta 



P 



P± 



P 



(14) 



in the rest frame P = (0, 0, 0, iM x ) of the pion and 
real p. To evaluate the quark DSE at these com- 
plex arguments it is convenient to determine the 
quark dressing functions A(p) and B(p) as func- 
tions of p instead of the usual p 2 . Numerically we 
solve the quark DSE for fixed gluon propagator and 
quark-gluon vertex by a plain fixed point iteration 
method. This iteration has to be carried out on 
the grid of real parts pj of p^ with fixed imaginary 
parts parametrised by Af T . Thus for any given 
value of M„ the quark DSE can be solved self- 
consistcntly without requiring knowledge of the 
propagator at a different M w . Of course, due to 
symmetries it is not necessary to iterate the DSE 
for the full grid of momenta p^ but for a subset 
parametrised by the absolute value of the pj , their 
four components and the pion momentum 

When solving the quark DSE, FigJTJ in the com- 
plex plane one encounters the following ambiguity: 
since analytic continuation is not unique anymore 
on the discretiscd Matsubara frequencies one ob- 
tains results in the complex plane which are de- 
pendent on the momentum routing in the dress- 
ing loop of the DSE. This ambiguity is dependent 
on the volume. It is small for volumes larger than 
V = (2fm) 4 , i.e. in the region where volume effects 
are small anyhow. Below V = (1.5 fm) 4 , however, 
there are sizeable effects. In this work we choose a 
momentum routing where the complex momentum 
goes through the quark part of the loop, whereas 
the gluon part depends on the real loop momentum 
only. Other choices probe the gluon propagator in 
the complex plane where no numerical results are 
available so far. 



D. Truncation scheme for the quark-gluon 
interaction 



We now proceed by specifying the explicit in- 
teraction used in the quark-DSE Eq. © and the 
BSE Eq. (TTU)) . It consists of two distinct parts, the 
dressed gluon propagator and the dressed quark- 
gluon vertex. For the gluon propagator we use 
solutions from the ghost and gluon DSEs as ob- 
tained in Refs. [22|, [23J. There the system of DSEs 
has been solved using a specific truncation scheme 
for the ghost-gluon vertex and the three-gluon ver- 
tex as input. The resulting dressing functions for 
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P [GeV] 

FIG. 3: Finite volume effects in thegluon dressing 
function Z(p 2 ) as determined in Ref. [22( [. 




FIG. 4: The approximated Schwinger-Dyson equation 
for the quark propagator with effective one-gluon ex- 
change and one-pion exchange. 

= if + + 1 

FIG. 5: The Bethe-Salpeter equation corresponding to 
the quark self-energy of Fig. [4] The up-down arrow 
indicates an averaging procedure of the pion exchange 
diagram with its counterpart where the upper pion- 
quark vertex is bare and the lower one dressed. 



the ghost and gluon propagator arc in good quali- 
tative agreement with corresponding lattice results 
with a quantitative difference of the order of ten 
percent at the peak of the gluon dressing function 
for momenta around p = 1 GeV. This is in marked 
contrast to corresponding dressing functions ob- 
tained from other approaches as e.g. the back- 
ground gauge formalism [3l| , where there is almost 
an order of magnitude difference between the char- 
acteristic scale of the DSE-solution and the lattice 
result. Thus, the results of [3l[ cannot be used in 
phcnomenological calculations as they would ne- 
cessitate an artificially strong quark-gluon vertex 
in order to provide for enough interaction strength 
to generate dynamical chiral symmetry breaking. 

In general, two types of solutions for the ghost 
and gluon propagators have been found in the deep 
infrared momentum region. There is a 'scaling' so- 
lution corresponding to an infrared vanishing gluon 
propagator and an infrared diverging ghost dress- 
ing function and there is a continuous family of 
'decoupling' solutions corresponding to an infrared 
finite gluon propagator and an infrared finite ghost 
dressing function, see e.g. H^, [32], HH and refer- 
ences therein. Recently, this family of decoupling 
solutions plus the limiting scaling one has been 
connected with ambiguities of fixing Landau gauge 
completely [34j |. It is furthermore a current issue 
of intense debate whether and how these two so- 
lutions relate to important fundamental questions 
like the mechanism of confinement of QCD and the 
possibility of a nonperturbative BRST-symmetry. 
Fortunately for the purposes of this work the be- 
haviour in the deep infrared is irrelevant since scal- 
ing vs. decoupling distinguishes the ghost and 
gluon at momentum scales below p 2 = 50 MeV 2 
corresponding to volumes far larger than the ones 
investigated here. 

For the gluon propagator at finite volume and 
discretisation we could in principle use results 



obtained on a torus [22|, cf. Fig. [3J How- 
ever, since corresponding results for the finite vol- 
ume quark-gluon vertex are not yet available we 
believe it could be more systematic if we ne- 
glect volume effects in the gluon propagator as 
well. The numerical results in sections IIII Al and 
IIII Bl are therefore obtained with the infinite vol- 
ume/continuum gluon propagator evaluated at the 
bosonic Matsubara frequencies appearing in the 
quark DSE Eq. ([5]). Nevertheless we also per- 
formed some calculations including finite volume 
effects in the gluon propagator which are discussed 
in section IIII B[ 

Next we give an explicit expression for the 
quark-gluon vertex. Here we follow the strategy of 
Ref. [19( and split the vertex into a part contain- 
ing only effects from pure Yang-Mills theory and 
a part containing the effects of pion backreactions 
to the quark propagator. The resulting interaction 
is given diagrammatically in figures 2] and [5] It 
satisfies the axial-vector Ward-Takahashi identity 
which is mandatory to obtain the pion as a Gold- 
stone boson in the chiral and infinite volume limit. 

The Yang-Mills part of the quark gluon vertex 
in this approximation is given by fl9j 



Tl M (k*) = lv Z 2 /Z z T YM (k 2 ) 

k 2 \ -1/2-* 



iYM / .2 



(15) 
(16) 



di 



+ 



d 2 + k 2 d\ + (p - d 2 f A| CD + fc 2 



4tt 



A 2 



\ V V QCD / 



k 2 — A 2 

QCD 



-28 



with the gluon momentum fc 2 , the one- loop value 
S = 44j^z8AT f° r ^ ne anomalous dimension of the 



G 



rc-summcd ghost dressing function and = 0.2. 
We also use AL n = 0.52 GeV 2 similar to the scale 
obtained in Ref. [35J . The infrared exp onent k 
has been determined analytically in |36l . |37| and 
is given by k = (93 - VT^Ol) /98 ~ 0.595. While 
Z2 is determined selfconsistently in the quark-DSE 
we choose Z3 = 1 for the ghost renormalisation 
factor. Since the combination of a^, the gluon 
dressing function and the two factors of I/Z3 from 
the bare and dressed quark-gluon vertex together 
form a renormalisation group invariant this choice 
is possible without affecting multiplicative renor- 
malisability of the quark-DSE. The only remaining 
free parameters of our interaction are d\, d 2 and 
c?3. These have been determined in (l9| to repro- 
duce the physical pion mass and decay constant as 
well as a reasonable mass of the 77' in the chiral 
limit [lU. The resulting choice is di = 1.45 GeV 2 , 
d 2 = 0.1 GeV 2 and d 3 = 3.95 GeV 2 . 

The construction Eq. (fT5]) follows the frequently 
used rainbow-ladder approximation of the full ver- 
tex and therefore involves only the 7^-part of its 
tensor structure. Note, however, that it has been 
shown in Ref. [38| that such a model cannot cap- 
ture all essentials of dynamical chiral symmetry 
breaking. Nevertheless it represents a useful start- 
ing point for our investigation. 

It remains to specify the quark meson vertex for 
the backreaction of the pion to the quarks. Here 
we employ only the leading part of the pion ampli- 
tude Eq. (fTTT) given by the function F 1 (p i ;P). In 
Ref. [19| the chiral limit relation Eq. (fT2]) has been 
used to represent the pion. However, since B = 
in the chiral limit on a torus this is not an option 
here. Instead we use the approximation 

r*(p t ;P) = 7 5^i(3?(p 2 ),K • P = 0,P 2 ) (17) 

where 5i(p 2 ) is the real value of p 2 . Since F\ is 
calculated from the volume dependent pion BSE 
this approximation takes volume effects in the pion 
wave function into account. As we will see, this 
model is elaborate enough to obtain a qualitative 
picture of volume effects due to the pion backreac- 
tion. 



III. NUMERICAL RESULTS 

Having discussed the details of our truncation 
for the quark-gluon interaction we now proceed to 
present our numerical results. All results shown 
are obtained using antiperiodic boundary condi- 
tions in all four space-time directions. We will also 
use L = T\ other choices are possible in principle 
and will be dealt with elsewhere. 



A. Volume and quark mass dependence of 
the quark condensate 

Let us first discuss the volume and quark mass 
dependence of the quark condensate. With finite 
volume V = L 3 T this quantity can be extracted 
from the dressed quark propagator via 



AZ 2 N C ^ 



L 3 T 4- p 2 A 2 (p 2 )+B 2 {p 



(18) 



The result is a renormalisation point independent 
but cutoff dependent condensate. At finite bare 
quark masses, necessary in a finite volume, this 
quantity contains quadratic divergences in the con- 
tinuum limit of infinite ultraviolet cutoff. There- 
fore in principle one has to think about the neces- 
sity of subtractions. However, it has been argued 
in Refs. 0,0] that subtractions are not necessary in 
the finite volume scaling regime, i.e. in the region 
of volume and quark mass where the condensate is 
a function of the dimensionlcss variable fi — mVT. 
only. Here £ represents the chiral condensate in 
the infinite volume limit. In this region the ul- 
traviolet divergences of the form mA 2 and m 3 In A 
with cutoff A are suppressed by 1/V and 1/V 3 re- 
spectively. These corrections are non-universal. 

The leading behaviour of (qq)(fx), however, is 
universal in the scaling regime and can be calcu- 
lated from random matrix theory (RMT) [f|. For 
quarks in the fundamental representation of SU (3) 
gauge theory the result is given by 

= n[iN s +v{n)K Ns +v{y) 

+lN f +u+i{p) Kn { +v-i(p)\ (19) 

where the I n and K n are modified Bcsscl functions 
and v is the topological charge. 

In the following we determine the quark con- 
densate in this scaling regime from solutions of 
our quark DSE for appropriate volumes and quark 
masses. We work in quenched approximation, i.e. 
Nf = 0. Our result is displayed in the left dia- 
gram of Fig. |5] together with RMT scaling accord- 
ing to Eq. (fH?]) , We clearly do observe scaling of 
the condensate with the dimensionless variable \i. 
However, the quantitative details of our results do 
not agree with the scaling behaviour predicted by 
RMT. So what is going wrong ? The reason for 
this discrepancy can be explained from the lim- 
itations of our truncation scheme: of course one 
would expect to see RMT-scaling if and only if 
all Green's functions used in the quark-DSE would 
scale with \x. Unfortunately this is not the case 
for our volume and quark mass independent in- 
put, namely the gluon propagator and the dressed 
quark-gluon vertex. These quantities constitute a 




FIG. 6: Left diagram: The quark condensate (qq) normalised by its infinite volume chiral value E as a function 
of the dimensionless variable /i = mVE. Shown are mean field scaling and the result from random matrix theory 
compared to our results from the quark-DSE for several lattice sizes and box lengths. 

Right diagram: The quark condensate evaluated as a function of box length L for two different quark masses m 
corresponding to pion masses of M w = 100, 300 MeV. 




L [fm] L [fm] 



FIG. 7: Left diagram: The pion mass as a function of box length L for two different infinite volume pion masses 
once without ('without tt') and once including the pion backreaction to the quark propagator. 
Right diagram: Same as left diagram but for the pion decay constant. 



/i-indcpcndcnt background in which the quark con- 
densate is evaluated. Nevertheless the volume and 
quark mass effects associated with the quark-DSE 
Eq. © itself are present. As a result one should 
expect mean field scaling of the quark condensate 
as a function of \i: 

(qq)(p)~H 1/S , (20) 

with 5 = 3. This is indeed what we observe. As can 
be seen from Fig. [6j in the scaling region fj, < 1 our 
result satisfies Eq. (|20l) with very good accuracy. 
We also observe critical slowing down in our nu- 
merics. We have checked that the scaling Eq. (|20|) 
persists in our solutions regardless of the details 
of the truncation of the quark-gluon vertex and 
also for different boundary conditions of the quark 
field. Around /i ~ 1 the aforementioned cut-off ef- 
fects set in and take over completely for /i > 10. 
It will certainly be interesting to check whether 



RMT-scaling appears in our approach when the re- 
strictions on the quark-gluon vertex are lifted and 
volume and quark mass effects in the vertex are 
taken into account. This is, however, beyond the 
scope of the present work. 



In the right hand diagram of Fig. [S] we also show 
the quark condensate evaluated as a function of 
box length for two different quark masses. Again 
we clearly observe that there is a critical region of 
volume below which the pattern of chiral symme- 
try breaking is changed substantially and finally 
lost for very small volumes [13, HH . This happens 
roughly around L — 1.5 fm which is also the char- 
acteristic length scale for the mean field scaling we 
observe. 
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FIG. 8: The RlM-n] ratio as a function of box length. Shown are our results for two different pion masses without 
(left diagram) and with the inclusion of pion cloud effects (right diagram) . The latter results are compared to 
corresponding ones from chiral perturbation theory The inlays magnify the region where our results have a 
negative sign. 



B. Volume and quark mass dependence of 
the pion mass and decay constant 

In the previous section we discussed the vol- 
ume dependence of a gauge invariant quantity, the 
quark condensate. Next we focus on two other 
gauge invariant quantities, the pion mass and the 
pion decay constant. In Fig. [7] we show our results 
for the volume dependence of these quantities. In 
order to also discuss effects due to different bare 
quark masses we include results corresponding to 
two different pion masses in the infinite volume 
limit, M„ w fOO, 300 MeV. Furthermore in each 
diagram we show results with and without taking 
the pion backreaction to the quark propagator into 
account. 

As concerns the pion mass, we see only very 
small effects for volumes larger than V = (1.8fm) 4 . 
These volumes seem to be large enough to be 
close to the infinite volume limit. For the pion 
decay constant, however, the situation is differ- 
ent: there are deviations from the infinite volume 
limit of the order of five percent up to volumes of 
V = (2.3fm) 4 . Recalling Eq. (13]) wc find that the 
calculation of /„• directly includes the pion wave 
function. On the other hand this quantity is absent 
in the rainbow-ladder calculation of the pion mass 
and only present at a sublcading level when we ad- 
ditionally backcouple the pions to the quarks. Our 
result for fir(L) then indicates that volume effects 
in the wave function of the pion are larger than the 
ones in its pole mass. This is as expected from the 
general considerations of Ref. [l| . 

In Fig. [S] we present results for the quantity 



R\MJ = 



M„{L)-M„{oo) 



(21) 



which is sensitive to the deviation of the pion 



mass in a box of length L compared to the infi- 
nite volume limit. In the left diagram we show 
results without the inclusion of the pion backreac- 
tion on the quark propagator, in the right diagram 
we show results taking these pion cloud effects 
into account. We study two different pion masses 
M n = 100, 300 MeV. In general the size of the vol- 
ume effects are larger for the smaller pion mass as 
expected. In both calculations we observe a posi- 
tive shift of i?[Af 7r ] if the box length is very small, 
followed by a turnover and a negative shift for box 
lengths larger than L ~ 1.6— 1.9 fm. The negative 
shift is larger and sets in at smaller volumes if pion 
cloud effects are omitted (left diagram) . Therefore 
it seems as if this negative shift may predominantly 
be a property of the quenched theory. Indeed, the 
very same sign change in the i?[Af^]-function has 
also been observed in quenched lattice QCD @. 
Our results without pion backreaction are there- 
fore in qualitative agreement with the lattice re- 
sults. In order to compare also quantitatively we 
would have to introduce periodic boundary condi- 
tions in the spatial directions as well as different 
box lengths in the space and time directions. This 
study is relegated to future work. 

In the right diagram of Fig. [3] we compare our 
results (including pion cloud effects) with a corre- 
sponding one for the heavy pion extracted from a 
rcsummed form of Liischcr's formula together with 
input from chiral perturbation theory at next-to- 
next-to-leading order [li|. For the light pion the 
xPT-calculation is outside the range of applica- 
bility of the p- counting- regime and therefore not 
shown in the plot. For the heavy pion mass our 
results show the onset of volume effects at similar 
length scales than \ however the sign is dif- 
ferent. xPT predicts a positive shift in R[Mtt] for 
all volumes. Since we have seen that our nega- 
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L [fmj 

FIG. 9: The chiral condensate as a function of box 
length once calculated with and once without finite 
volume effects in the gluon propagator. 

tive trend in i?[M„.] indeed gets smaller when we 
take the pion backreaction into account we can not 
exclude that our result even changes into the pos- 
itive domain when we treat the pion cloud effects 
more systematically. This will be explored in fu- 
ture work. We wish to note, however, that nega- 
tive values for i?[M w ] have also been observed in 
the quark meson model for unequal spatial and 
temporal box lengths, see [TtJ for details. 



C. Additional volume effects due the gluon 

Finally we discuss additional volume effects in 
the chiral condensate due to volume effects in the 
gluon propagator. As discussed in section IIIDI 
these effects should be included together with the 
corresponding effects in the quark-gluon vertex. 
Nevertheless it may be instructive to explore the 
size of the effects due to the gluon alone. For the 
quark mass function these have already been dis- 
cussed in Rcf. dHH]. Note that in [H[ the volume 
effects in the gluon propagator have been overes- 
timated and consequently unrealistically large ef- 
fects in the quark mass have been found. This find- 
ing has been corrected in [25| using the updated 
calculation of the finite volume gluon propagator 
from [22j . Here we use a similar input, as detailed 
in section 111 Dl 

Our results are shown in Fig. |H1 We clearly 
observe a similar qualitative behaviour with and 
without the finite volume gluon. The quantitative 
differences between the two calculations are small 
(i.e. below five percent) apart from a small region 
around L = 1.5 fm, where the vanishing of chiral 
effects happens at slightly different scales. For very 
small volumes both results become equal. This is 
easily understood since in this region the quark 
condensate is dominated from the bare quark mass 



and dynamical effects are almost absent. We also 
performed calculations of the pion mass and decay 
constant using the finite volume gluon propagator 
and found even smaller effects there. Therefore we 
agree with the previous finding of [25[ : the effect of 
finite volume corrections in the gluon propagator 
on physical obscrvables is small. 



IV. SUMMARY AND CONCLUSIONS 

fn this work we presented an exploratory study 
of volume effects in the quark condensate, the pion 
mass and its decay constant determined from a 
coupled set of Dyson-Schwinger (DSEs) and Bcthe- 
Salpeter equations (BSEs). These equations are 
capable of connecting the large and small regions 
of momenta, quark masses and volumes, and there- 
fore provide an interesting tool to determine these 
effects systematically. In our exploratory study we 
concentrated on volume effects generated in the 
quark DSE and the pion BSE and neglected cor- 
responding effects in the quark-gluon vertex. This 
deficiency introduces a volume independent scale. 
Consequently, at small volumes we do not observe 
the universal scaling of the condensate func- 
tion of fi = mVE predicted from random ma- 
trix theory. Instead we find critical scaling with 
mean field exponents at box lengths of L « 1.5 fm. 
In general we observe strong volume effects in all 
quantities for box lengths below L = 1.8 fm. For 
the pion these effects are somewhat larger in the 
wave function than in its pole mass. As a conse- 
quence the effects in the pion decay constant 
are also larger than in the pion mass. In the re- 
gion above L = 2 fm we find volume effects with 
opposite sign than the ones determined by chiral 
perturbation theory. Whereas chiral perturbation 
theory predicts an increase of the pion mass we 
find a slight decrease for L = (2 — 4) fm with the 
increase setting in only below L = 2 fm. In the 
quenched theory this behaviour is in qualitative 
agreement with results from lattice QCD 0. For 
the unquenched theory, however, the persistence of 
this effect in contrast to chiral perturbation the- 
ory points towards deficiencies in our truncation 
scheme. We expect this problem to be cured by 
taking into account explicit volume effects in the 
quark-gluon vertex along the lines of recent work 
on the infinite volume case [3(| [H| . This will be 
the subject of future work. 
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